Casimir Energies of Cylinders: Universal Function 

E. K. Abalc! and K. A. Miltorfl 

Homer L. Dodge Department of Physics and Astronomy, 
University of Oklahoma, Norman, OK 73019 

L. KaplanS 

Department of Physics, Tulane University, New Orleans, LA 70118 

(Dated: August 30, 2010) 

Abstract 

New exact results are given for the interior Casimir energies of infinitely long waveguides of 
triangular cross section (equilateral, hemiequilateral, and isosceles right triangles). Results for 
cylinders of rectangular cross section are rederived. In particular, results are obtained for interior 
modes belonging to Dirichlet and Neumann boundary conditions (TM and TE modes). These 
results are expressed in rapidly convergent series using the Chowla-Selberg formula, and in fact 
may be given in closed form, except for general rectangles. The energies are finite because only 
the first three heat-kernel coefficients can be nonzero for the case of polygonal boundaries. What 
appears to be a universal behavior of the Casimir energy as a function of the shape of the regular or 
quasi-regular cross-sectional figure is presented. Furthermore, numerical calculations for arbitrary 
right triangular cross sections suggest that the universal behavior may be extended to waveguides 
of general polygonal cross sections. The new exact and numerical results are compared with the 
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I. INTRODUCTION 

Casimir's seminal 1948 paper [lj, demonstrating that two perfectly conducting parallel 
plates placed in vacuum attract, created something of a stir in the physics community 
because of its unexpected nature. In principle, it was verified experimentally in 1958 by 
Sparnaay [2j, although conclusive measurements had to wait until the end of the 20th century 
3-5]. Even more sensational was Boyer's 1969 [6| calculation for the Casimir energy of 
a perfectly conducting spherical shell, which produced a repulsive force, rather than the 
attractive force that Casimir had postulated earlier 7| for a model of the electron. Since 
then, the question of how the sign and magnitude of the Casimir energy depend on the 
geometry of the system has been a constant leitmotif in Casimir energy studies. 

In this paper, we attempt to elucidate this question by calculating the Casimir energies for 
the interiors of infinitely long cylinders of various cross sections. Such calculations have long 
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11] cross sections, but 



been available for waveguides of rectangular SB and circular 
it apparently was not widely appreciated that energies for certain triangular cross sections 
were equally well calculable. We will work with triangular (equilateral, hemiequilateral, 
and isosceles right triangles), square, and rectangular cross sections. The triangular and 
rectangular cylinders are ideal since their eigenvalues are known explicitly. We obtain results 
in a very efficient form extremely amenable to numerical evaluation by use of the Chowla- 



Selberg formula 12Hl4|. which is far more rapidly convergent than the direct evaluation 



of Epstein zeta functions. Even better, we are able to obtain closed-form results for the 

n 

interior Casimir energies by use of a convenient representation found by many authors [15J- 
221 ] . These powerful methods supersede the previous method of directly summing the Epstein 
zeta function. With the obtained data, we can then relate Casimir energies for these cross 
sections with their distinguishing geometrical properties, notably the perimeter P and the 
cross-sectional area A. Our chief result is the discovery that there appears to be universal 
functional dependence for the Casimir energy per length multiplied by the area, SA, due 
to the interior fluctuations inside an infinitely long waveguide, in terms of a dimensionless 
quantity expressing the waveguide's cross-sectional attributes, A/P 2 . This universality class 
applies to figures based on regular polygons; for a rectangle, for example, of sides a and b a 
different dependence obtains, which changes sign, as expected, when the rectangle becomes 
sufficiently elongated. 



2 



It is interestin g to note that the square and circular waveguides have been extensively 



studied 
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3M27I|. in contrast to the triangular case which has only partially been ex- 
amined. As far as the authors can tell, this paper is the first with explicit numerical results 
for the interior Casimir energy for infinitely long waveguides of triangular cross sections. 
Earlier work has, however, been done on plane equilateral triangular geometries l28|. For 
other work on equilateral triangular domains see Refs. 
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3d. There has been recent 



work in which explicit eigenvalues for the equilateral triangle have been used in computing 
numerical values for a piston of such a cross section 31 ]. 

Probably the reason why these Casimir energies were never evaluated, although the ex- 
plicit eigenvalues appear in textbooks 32|, |33J, is that most authors considered interior 



calculations suspect. It is, of course, not possible by these methods to extract exterior 
eigenvalues or the corresponding Casimir energies. For smooth surfaces, such as a circular 
cylinder, only the sum of the interior and exterior Casimir energies of an infinitesimally thin 
boundary shell can be calculated unambiguously. In the language of heat kernels, the 
heat kernel coefficient is proportional to the cube of the curvature, and cancels only when 
both interior and exterior modes are included. When a2 ^ and only the interior is con- 
sidered, there is a logarithmic divergence that cannot be removed, leaving an ambiguous 
Casimir energy. But, for polygonal boundaries, there is no such coefficient, and the interior 
Casimir energies seem well defined. There are divergences associated with the corners, but 
they do not contribute to the global Casimir energy of the waveguide, because they may 
be unambiguously subtracted off. For a review of the state of the art for self-energies of 
smoothly bounded regions, see Ref. [34] 

The outline of this paper is as follows. In Sec. |TT]we formulate Casimir energies in terms 
of mode summations, which we regulate either by dimensional continuation or by a temporal 
point-splitting cutoff. The former approach leads naturally to the use of the Chowla-Selberg 
formula, while the latter produces an expression that is summed through use of the Pois- 
son summation formula; either approach can lead to a result that can, in special cases, be 
expressed in closed form. These machineries are applied to cylinders with equilateral trian- 
gular cross sections in Sec. IHIt for both Dirichlet and Neumann scalar modes, and hence 
for interior electromagnetic modes for a perfectly conducting boundary. The same is done 
in Sec. IIVI for a bisected equilateral triangular cross section, that is, a hemiequilateral or 
30°-60°-90° triangle. A square cylinder is reconsidered in Sec. |Vl as is a rectangular cylinder 
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in Sec. IVII The bisection of a square waveguide, one formed by a right isosceles triangle, 
is treated in Sec. IVII I General right triangular cross sections are considered numerically in 
Sec. IVHH showing that a finite Casimir energy is obtained even when the eigenvalues are not 
explicitly available. All the energies for triangular Dirichlet boundaries found in this paper 
lie on a universal curve, which smoothly joins the proximity force approximation. Why the 
interior modes of a continuous curved cross section, such as that of a circle, cannot yield a 
finite Casimir energy is discussed in Sec. IIXI The corner divergences are reflected only in 
the a\ heat-kernel coefficient, unlike the curvature divergences which show up also in the 
a 2 coefficient. Again, unlike for continuously curved figures in a plane, Casimir energies for 
squares and equilateral triangles, and figures obtained from these by bisection, are likewise 
unambiguously calculable, and given in Sec. [Xj along with numerical results for right trian- 
gles. Results for the cylindrical geometries considered are summarized in the Conclusions in 
Sec. IXII The Appendix sketches the derivation of the Chowla-Selberg summation formula 
from the Abel-Plana formula. 

II. CASIMIR ENERGY OF A CYLINDER 
A. Mode summation 

The Casimir energy per length, £, for infinite cylindrical geometries can generally be 
expressed as a sum over mode frequencies in the following form, 



where 7^ n are the eigenvalues of the two-dimensional Laplacian, (V^ + 7^ n ) $ mn (r±) = 0, 
subject to $ m „(r_i_) satisfying appropriate boundary conditions on the cylinder's surface. 
Here, r_|_ denotes the two coordinates transverse to the longitudinal cylinder axis. For 
cases where the eigenvalues have explicit formulas, 7^ n is a quadratic function of the mode 
numbers m and n. 

The form of this quadratic function depends on the geometry, and the range of the mode 
numbers depends on the boundary conditions being applied. In the subsequent sections, we 
will look at the cases of Dirichlet and Neumann boundary conditions. The electromagnetic 
case can be obtained from these two classes of modes. Because of the cylindrical geometry, 




(2.1) 
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the modes of the electromagnetic field in a perfectly conducting waveguide can be expressed 
in terms of two scalar fields satisfying separately Dirichlet (E or TM modes) and Neumann 



boundary conditions (H or TE modes) 32j, |33J. It follows then that the Casimir energy for 



an electromagnetic field is simply the sum of the Casimir energies for the E and H modes. 
(This breaks down if the boundaries are imperfect, for example, for a general dielectric 
cylinder [351].) 

The expression above for the energy is formally divergent and as usual requires a regular- 
ization technique in order to obtain physically interpretable results. Out of the many ways 
to regularize such an expression, we choose the methods of dimensional regularization and 
point-splitting regularization. 

B. Regularization methods 

1. Dimensional regularization 

By extending the one dimension for the longitudinal wavevector k to d dimensions, dk — > 
d d k, we obtain the following expression for the energy per length, 

r(-(l + d)/2) x - , 2 x(i+d)/2 , v 

g = -j™ 2 d +27r(rf+ l)/ 2 l^Umn) ■ (2-2) 

m,n 

In the limit d — > 1, we obtain a r(— 1) divergence. However, one can surmount this problem 
by analytic continuation. For the cases considered in this paper, the eigenvalue expression 
7^ n is a simple quadratic form in m and n. This allows us to employ the Chowla-Selberg for- 



mula, an exact formula for a class of Epstein zeta functions 12|, ll3j, proved in the Appendix, 



and thereby utilize the reflection property of the zeta function, 

r(a)C(2a) = r((l - 2a)/2)C(l - 2s)^- 1 l\ (2.3) 

which allows us to continue to d — > 1. The technique of dimensional regularization is our 
method of choice throughout this paper since the Chowla-Selberg formula converges expo- 
nentially faster than the Epstein zeta functions obtained by the point-splitting regularization 
we present below. 

In fact, it has been known for many years that we can sum the zeta function into a 



form 
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21] that can yield a closed-form expression for the Casimir energy in the cases of 



a square, isosceles right triangle, equilateral right triangle, and a 30°-60°-90° triangle. 



2. Point- splitting regularization (Cutoff) 



A more physical approach is to start from the general form for the quantum vacuum 
energy in terms of the Green's function for a given frequency uj, 

^s/w/s 2 "*^ (24) 

which for the situation of a cylindrical waveguide, gives the energy per length 

where in the second line we have made a Euclidean rotation, u — > i£, as well as introduced 
a time-splitting regulator, t — t' — >■ rr, where we are to take r to zero at the end of the 
calculation. Assuming r > 0, we do the ( integration by closing it in the upper half plane, 
with the result 

m,n J 00 

= I l im f-±) — r e -^. (2.6) 
2-oV drJJ_ 00 2n^ n 

Just as the dimensional regularization method was paired with the analytic continuation 
of the exponent, followed by use of the Chowla-Selberg formula, the cutoff regularization 
method is combined with Poisson's sum formula. After rearranging the summations to span 
—00 to +00, the summation is converted by Fourier transformation from one involving pow- 
ers of e _r to powers of e _1 / T . The final expression after simplifications is similar to the 
Chowla-Selberg result except for a slowly converging double sum. As noted above, in fact, 
that double sum can in special cases be summed into an explicit closed form; both regular- 
izations give the same answer and thus serve as somewhat independent consistency checks. 
In addition, the cutoff regularization, unlike the dimensional-zeta-function regularization, 
exhibits the explicit divergences. We are then able to identify a "volume" and a "surface" 
divergence, as well as a divergent term independent of the scale of the cross section. All of 
these divergent terms can be regarded as contact terms, and, therefore, may be discarded. 
Corner divergences occur only in the constant term. 
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a 

FIG. 1: Equilateral triangle of side a and height h = 



III. EQUILATERAL CYLINDER 

We consider an infinitely long cylinder of equilateral triangular cross section with side 
length a and height h, as shown in Fig. [TJ The eigenvalues of the Laplacian on an equilateral 



triangular surface with Dirichlet or Neumann boundary conditions (solved by G. Lame 36]) 



are of the form 7^ nj) = |p- ( m2 + n 2 +p 2 ), where m + n+p = 
the previous expression can be simplified resulting in 

4vr 2 



32 



331] . Using this constraint, 



TL= ^{m l + mn + n l ). (3.1) 

Doundary 



33]) have 



Here, m and n are positive or negative integers; their allowed values depend on the 
conditions, as indicated below. The eigenf unctions (which are given in Refs. 32L 
degeneracy under cyclic or anticylic permutations of the three indices m, n, p, so a factor 
of 1/6 must be inserted before the mode sum. One would, indeed, expect only two mode 
numbers as it is a two-dimensional surface; however, the earlier form is more useful for 
counting degeneracies. 



A. Dirichlet boundary 

Dirichlet boundary conditions require the fields to vanish on the boundary C, &{r±)\c = 
0. The mode numbers are now restricted to m ^ 0, n / 0, and p ^ 0. If any index 
vanishes, the corresponding eigenfunction would likewise vanish. One can then use Eq. (12. ip 
with Eq. (13. ip . and perform the appropriate sum over m and n to obtain the Casimir energy. 
The dimensionally regulated expression ( 12. 2 p may be evaluated by using the Chowla-Selberg 
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formula, for Res > 1: 



£ "{am 2 + 6mn + cn 2 y s = 2a~ s C(2s) + ^^7/^ (2^ - l)r(a - 1/2) 



m,n=— oo 



+ r(a)A a /2-V4 /^ 2Z Tia " 1/2q 'i-2 a Wcos(w7r6/o).R:i / 2- 8 (7i7rv / A/a). (3.2) 

^ ' * n=l 

However, it is valid only when the discriminant A = 4ac — b 2 > 0. The double prime 
indicates m = n = is excluded from the summation range. The divisor function, cjfc(n), is 
the sum of the k-th powers of the divisors of n, 

a k (n) = J^d"- ( 3 - 3 ) 

d\n 

In the particular case of the Dirichlet equilateral triangle, A = 3. The special mode 
sum to be evaluated involves the explicit exclusion of the case when one of the mode labels 
vanishes: 

r(-(l + d)/2) VA '- 2l(d+1)/2 



4 /7T\2 

3 \h 



Eq ' 2 a! + 2 7r( d+1 )/ 2 

xi ij2"(^ 2 + n 2 + mnf +d ^ 2 -6((-l-d) }. (3.4) 

V m,n 



After a few manipulations such as use of the reflection property (12. 3p we obtain 

4? = IT^( 8 <(3) -3 3/2 C(4) - 4vr 2 (12) 3 / 4 f;n- 3 / 2 (-l)V3(n)K3 /2 (n 7 rv / 3)). 



n=l 



(3.5) 



The above expression converges very fast and reaches an accuracy of twenty decimal places 
when just the first eight terms in the n sum are included, 

MD) 0.017789138469130117062 

E k = ^2 • (3- 6 ) 

B. Neumann boundary 

For Neumann boundary conditions, the normal derivative of the eigenfunctions must 
vanish on the boundary, <9 n $(rjJ|(7 = 0. The restriction on the mode numbers is less 
severe than the Dirichlet case: here one of the mode numbers can be zero, but not two 
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simultaneously. This gives a different summation range from that for the Dirichlet case, and 
after similar regularization manipulations, we obtain the result 

F (N) _ f (D) _ C(3) , , 

^Eq - ^Eq ^2 ^' a ) 

0-045982 

= ^— . (3.7b) 

C. EM perfectly conducting boundary 

As discussed in the Introduction, the electromagnetic Casimir energy in the interior of 
perfectly conducting waveguide having an equilateral triangular cross section is simply the 
sum of the two previous energies, 

47 =C+4: ) (3.8a) 
°' 028193 (3.8b) 



h 2 

Note that while the Dirichlet energy is positive, the electromagnetic energy is negative, 
because the H mode overwhelms the E mode. 



D. Point-splitting regularization 

Poisson's summation formula gives an expression equivalent to the original sum in terms 
of the summand's Fourier transform: 



E f(P)= E / dpe^fip)). 



(3.9) 



p=— oo q=—oo 

This is useful for converting formulas most readily applicable for a large parameter to one 
most convenient when the parameter is small — that is, it is a duality transformation. 

We will demonstrate this regularization method with the Dirichlet equilateral triangle. 
The mode sum in the energy (12.61) can be given explicit form as 



m,n \m,n=—oo m J 



(3.10) 



where 7^ n is given by Eq. (I3.ip . and the sums on the right hand side extend over all integer 
values, positive, negative, and zero. We first perform a transformation on 7 mra with m = r + s 
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and n = r — s. It follows that 



m,n 



+ e 



r+l/2,s+l/2 



3 ^( e -V fc2 +7o 2 ,) + 2e- rV ^V 



(3.11) 



where now 7 2 s = |p-(3r 2 + s 2 ), and the r and s sums again extend over all integers. In a 
manner similar to Ref. js], we adopt spherical coordinates and evaluate the generic Fourier 
transforms: 



dk Yl e- T ^ k2+ P {a{r+b)2+c{s+d)2} 

r,s 



ibp —2mdq 



p,q 



OO /•■7T 

2 



-4= / iftiH / sin0d0 
PVacJo Jo jo 



2n 



(3.12) 



where we have defined the vectors R = (^/J3ar, y/JJcs, k) and p = -^{p/y/a,q/y/c, 0). Car- 
rying out the radial and angular integrals, we find straightforwardly 



dk 



k 2 +/3{a(r+b) 2 +c(s+d) 2 } 



r. s 



E 



p.q 



-2m 



i(bp+dq)p 



iTYT 



/oo 



ac [r 2 /3 + 47r 2 (p 2 /a + g 2 /c)] 2 ' 

(3.13a) 



/3 



2ttt 



/ V« (r 2 /3 + 47rV/a) 



2/^3/2' 



dk e 



2 
r 



(3.13b) 
(3.13c) 



in terms of the abbreviation — | (f) 2 - With the appropriate coefficients for the Dirichlet 
equilateral triangle case, we isolate the divergences in the energy: 



= Km f - + J_ 

Eq r^O 27T 2 T 4 4ttt 3 6ttt 2 



lim 



3A 



P C 
+ 



o V 2vr 2 r 4 8vrr 3 48vrr 2 / ' 



(3.14) 



We note that the "volume" and "surface" divergent terms, which are respectively propor- 
tional to the area of the triangle A = h 2 / \/3 and the perimeter P = 2y/3h, are as expected, 
and are presumably not of physical relevance. The last term, a constant in h, certainly 
does not contribute to the self-stress on the cylinder. Only this term reflects the corner 
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divergences, and for the equilateral triangle we have C = 8. For a general polygon, with 
interior angles a,, the corner coefficient is 



7T a,; 

The coefficients in Eq. (I3.14p are proportional to the heat kernel coefficients discussed in 
Sec. HXl below. 

After these terms are isolated, the finite part of the energy is given by the following 
expression: 

^ = 144^ [ 12<{3) ~ ~ C(4) " — L i (p 2 + g 2 /3)2 j • ( 3 - 16 ) 

The double sum in the above expression converges slowly. Summing the first 500 terms 
for m and n, we reach a seven- decimal accuracy, 

C - (3.17) 

which agrees with Eq. (13. 6p . Because of the numerical agreement, one is convinced that 
both regularization methods yield exactly the same answer and we have thus shown that 
our calculations are correct. 

A virtue of the alternative form (13.161) is that it is now evident from Eq. (13 . 7aj) that 
in the electromagnetic energy (13.8aj) the ((3) term completely cancels. As a result, the 
electromagnetic energy is manifestly negative. This is a feature that will persist in all the 
examples treated in this paper. 



E. Closed-form result 



Remarkably, for the integrable polygonal figures we are considering, the Casimir energy 



can be given in closed form. Following Refs. 
cross section, from Eq. ( 12.21) . 



17H21|. we write for the equilateral triangular 



CD) 
Eq 



ton i(4rr)T( a ) ) [C(*)LsOO " C(2a)) 



in terms of the single series 



I*® = J2 



n=0 



1 



(3n + l) s (3n + 2) 



(3.18) 



(3.19) 
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a/2 



FIG. 2: 30°-60°-90° triangle obtained by bisecting an equilateral triangle. 



Now, this function satisfies the reflection property 

L a (s )V[s) = -L I — s). 

iy 1 v 1 2sins7r/2 dV ; 

Then, using Eq. ( 12 . 31) . we can take the limit s — > — 1: 

1 



(3.20) 



96/i 2 



^(2) - -C(3) 



7T 



(3.21) 



In fact, the remaining sum has a closed form: 



L 3 (2) = - [V'(l/3) - ^(2/3)] 



(3.22) 



in terms of the polygamma function. Thus 

1 



96h 2 



^ [^(1/3) - ^(2/3)] - -C(3) 



0.0177891 



(3.23) 



It is a priori remarkable that such an explicit form can be achieved for a strong-coupling 
problem. 

In particular, from Eqs. (13.7aj) and (I3.8al) . we see that the interior Casimir energy for a 
perfectly conducting cylinder with equilateral triangular cross section has the simple form 



c(EM) 

Eq ~432/^ 
that is, as we noted above, the £(3) cancels. 



[^(1/3) -V'(2/3)], 



(3.24) 



IV. HEMIEQUILATERAL CYLINDER 

The hemiequilateral triangle or 30°-60°-90° triangle is appropriately named since it is 
obtained by bisecting the equilateral triangle, as shown in Fig.[2j Hence, taking an equilateral 



12 



triangle of side length a and height h we obtain a hemiequilateral triangle of height h and 
side lengths a and a/2. This relation between the two triangles proves very useful in realizing 



that the hemiequilateral modes are a subset of the equilateral modes |32J. (The proof that 
there are no additional modes is given in that reference.) Therefore, one simply has to select 
the equilateral modes that satisfy chosen boundary conditions on the bisector. It is implied 
that the eigenvalues are of the same form as the equilateral triangle's, Eq. (13. ip . 



A. Dirichlet boundary 

To satisfy Dirichlet boundary conditions for the hemiequilateral triangle, the eigenmodes 
must vanish on all three sides of the hemiequilateral triangle. The Dirichlet modes for the 
equilateral triangle already satisfy that condition on two sides, so we are to select the ones 
that vanish on the bisector. After a few manipulations, the Casimir energy for the Dirichlet 
hemiequilateral triangle is related to the equilateral Dirichlet result with the formula, 

MD)_4° ) C(3) _ 0-0567229 

a positive energy again. 

B. Neumann boundary 

Similarly, the Neumann modes for the hemiequilateral triangle have to be chosen from the 
Neumann equilateral modes such that their normal derivative on the bisector vanishes. With 
this new restriction on the mode numbers, the relation between Neumann hemiequilateral 
and equilateral Casimir energies is established, 

AN) _ 4f _ C(3) _ 0.0708193 

t3m ~ 2 8vr/i 2 " ' 1 ' 

which is again negative. 



C. EM perfectly conducting boundary 

Straightforwardly, we obtain the electromagnetic Casimir energy in the interior of an 
infinitely long cylinder of hemiequilateral triangular cross section by adding the E and H 
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FIG. 3: Square of side a. 



modes, 



c(EM) _ AD) AN) _ 1 -(EM) _ 0-0140964 
°369 ~ °369 369 ~~ 2 q ~~ h 2 ' 



which is, remarkably, exactly one-half that of the energy of equilateral triangle, as might 
have been anticipated naively. 



V. SQUARE CYLINDER 



For completeness, we consider a square of side length a, as illustrated in Fig. El Of course, 
this is the geometry originally considered by Lukosz (sh and by Ambj0rn and Wolfram Jq]. 



The form for the eigenvalues is 7^ n = ^(m 2 + n 2 ) [32|, [33|. It is of a simpler quadratic 
form than the previous two geometries, which makes the counting of degeneracies somewhat 
simpler. 



A. Dirichlet boundary 

The Dirichlet modes for the square must vanish on all four sides, and so are products of 
two sine functions, a classic result in most electricity and magnetism books. The eigenvalues 
must therefore both be positive, m > and n > 0. With the appropriate summation range 
and dimensional regularization, we obtain a result consistent with that found in Refs. jsl, S|, 
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where we also display the result found by the alternative cut-off method: 



(D) 
Sq 



327r 2 a 2 
f 



2C(4) - vrC(3) + Svr 2 J2 r^ 2 a 3 (l)K 3/2 (nl) 



1=1 

oo 



327r 2 a 2 ^ 
0.00483155/a 2 



<( 4 )- 2 <( 3 )+ 4 E(]^p)2j 



(5.1a) 

(5.1b) 
(5.1c) 



The Chowla-Selberg formula ( 15. lap is extraordinarily convergent. 

Again this result can be given in closed form. Actually, we can do this directly from the 
double sum 



15 



16 



181 



where 



oo 

L *( 2 ) = E (O = G > 



where G = 0.915966 ... is Catalan's constant. Then 



£, 



(D) = 1 

Sq 167ra 2 



7T 



C(3) - -G 



0.00483155 



(5.2) 



(5.3) 



(5.4) 



B. Neumann boundary 

The Neumann modes for the square must have their normal derivative vanish on the 
boundary. Such modes are the products of two cosine functions, which implies that m and 
n cannot both be zero. The Neumann energy in relation to the Dirichlet energy is 



C(3) 



which again reproduces the known result. 



C. EM perfectly conducting boundaries 



0.0429968 
a z 



(5.5) 



Summing Dirichlet and Neumann energies, we obtain the following electromagnetic 
Casimir energy for the interior of a perfectly conducting square waveguide, expressed only 
in terms of one transcendental number G: 

G 0.0381653 



(EM) 
Sq 



AD) AN) 
C Sq C Sq 



24a 2 



(5.6) 
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FIG. 4: Isosceles right triangular waveguide, of side a. 



VI. RECTANGULAR CROSS SECTION 



The more general expression for a Dirichlet rectangle is easily obtained by the above 
methods. Again we display the result in the two forms, 

1/2 oo 



(D) 
Rect 



32vr 2 a6 
1 

32vr 2 a 2 



h f h\ * . 

-tt-C(3) + 2C(4) + Svr 2 - Yl l- 3/ Wl)K 3/ 2(27ilb/a) 
a 1=1 



(6.1a) 



1 + 



<(3) 



m,n=l 



m 



ax 2 
n 6 



-2 



(6.1b) 



For a = 6, a square, the Dirichlet energy is positive, but when one side is much larger than 
the other, the sign of the energy must change, for that situation corresponds to the classic 
case of Casimir attraction. Indeed, the Dirichlet Casimir energy vanishes for b/a = 1.74437, 
and is negative for larger values of b/a. See Fig. |9] below. Note that for a general rectangle, 
a 7^ b, a closed- form expression for the energy is apparently not achievable 19]. 
The corresponding Neumann result is simply 



(TV) 
Rect 



(D) 
Rect 



C(3) 



(6.2) 



167T0 2 

which again shows from Eq. (16.1bj) that the electromagnetic Casimir energy is always nega- 
tive. 



VII. ISOSCELES RIGHT CYLINDER 

The relation of the isosceles right triangle to the square is analogous to that of the 
hemiequilateral triangle to the equilateral triangle. Bisecting a square produces an isosceles 

16 



right triangle, as seen in Fig. HI The eigenmodes of the isosceles right triangle are a subset 
of the square eigenmodes satisfying the extra condition that they obey given boundary 
conditions on the corresponding square's diagonal. The eigenfrequencies are therefore of the 



same form as the square's, % mn = ^(m + n ), with different ranges for m and n 
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A. Dirichlet boundary 

Dirichlet boundary conditions on the isosceles right triangle translate into the constraint 
m > n, and n > 1. Just as in previous cases, the result follows seamlessly. The Casimir 
energy can be related to that of the Dirichlet square, 



,(D) 

a positive number again. 



C (D) g V C(3) 0.0263299 
tls ° ~ 2 + 16vra 2 ~ a 2 ' 



B. Neumann boundary 

By imposing Neumann boundary conditions, the mode numbers must now satisfy m > n, 
and n > 0. The relation to the Neumann case for the square is then, 



,(JV) 

again a negative number 



AN) _ 4/ _ m_ 0.0454125 

^ Iso ~ 2 16vra 2 a 2 ' 1 J 



C. EM perfectly conducting boundary 

We combine both results to obtain the electromagnetic Casimir energy for a perfectly 

conducting waveguide of right isosceles triangular cross section, 

~(em) _ AD) AN) _ lp(EM) _ G 0.0190826 

^iso — c iso +^iso — 2 s q — 48 a 2 — a 2 ' v'-Jj 

again a remarkably simple and unexpected result. 



VIII. GENERAL RIGHT TRIANGULAR CYLINDERS 



As a generalization of the hemiequilateral triangular cylinder and isosceles right triangular 
cylinder, we may consider cylinders having as their cross sections arbitrary right triangles 
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with hypotenuse a and angles 9, tt /2 — 9, and tt/2. Of course, away from the special 
values 9 = ir/6, 7r/4, or ir/3, it is not possible to obtain explicit expressions for the triangle 
eigenfrequencies 7 n . Instead, we may employ the efficient scaling method introduced by 
Vergini and Saraceno {stJ and further developed by Barnett jssj to obtain numerically all 
Dirichlet eigenfrequencies 7 n for a given triangle with high accuracy, up to some desired 
maximum frequency 7 max . The sum over eigenfrequencies is then evaluated numerically 
using point-splitting regularization as in Eq. (12. 6p . 



£,Ur) = U° 0d £ E V^T^e-v^, (8.1) 



where the error associated with terminating the sum at 7 max is proportional to e~ T7max , and 
thus we must consider T7 max 3> 1. 

To extract a physical Casimir energy at small r, we must first subtract off the known 
divergences associated with the first three terms in the Weyl expansion for the spectral 
density, associated with the area, perimeter, and corners of the triangular cross section, 

™ - s /: s r @ - a + ^ <-> 

where for the right triangle A = ^ sin(2#), P = a(l + cos # + sin #), and the corner coefficient 
is given by Eq. (13. 15ft . C = ^^J^ ■ For 7 max — > oo at fixed r, the divergent term recovers 
the form given earlier in Eq. (13.141) . 

We notice that even after subtracting the divergent terms, the limit r — > may not 
be taken numerically due to the constraint T7 max ^> 1, imposed by the fact that we have 
included only a finite part of the eigenfrequency spectrum. Instead we fix 7 max and fit 

£num(r) - £ d iv(r) = £ + ZiT + z 2 t 2 In T + z 3 t 2 H , (8.3) 

over a range r 6 [ti,t 2 ], where 7~^ x <C T\ < t 2 a, and where on the right hand side 
of Eq. (18. 3 p we have listed all non-divergent terms that may appear for a generic cavity, 
through order r 2 . The finite term £ is the Casimir energy per unit length for the cylindrical 
waveguide of chosen cross section. 

For the numerical results shown in Fig. El we have set a = 1, 7 max = 400, t\ = 0.05, 
and T2 = 0.2. With these parameters, the numerical approach reproduces exact results for 
the hemiequilateral and right isosceles triangles to better than 1%. Clearly, the numerical 
method briefly discussed here is not restricted in any way to right triangular cylinders or to 
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FIG. 5: (Color online) Graph of the interior Dirichlet Casimir energy per length E for various 
cylinders, multiplied by the area A. These are plotted as a function of the dimensionless ratio 
of the area A to the square of the perimeter P, A/P 2 . The points shown are for the 30°-60°- 
90° triangle, the isosceles right triangle, the equilateral triangle, the square, and the circle (which 
includes interior and exterior contributions), as well as other right triangles evaluated numerically. 
The curve is the PFA approximation, which should become exact when one of the angles of the 
triangle goes to zero. 

Dirichlet boundary conditions, but may be applied to arbitrary cavities in any dimension, 
provided that a reliable method of numerically evaluating the spectrum is available. A very 
similar approach was used in Ref. [39( to evaluate the vacuum energy in quantum graphs 
where no explicit solution is available. 

The graph shows the comparison of the exact and numerical data to the proximity force 
approximation (PFA) for the triangles, in which the energy is computed as though each 
opposite element of a right triangle, with the acute angle 9 — > 0, has the interaction energy 
given by that for parallel plates. This gives the formula 





This result closely matches the numerical data for very acute triangles. 
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IX. CIRCULAR CYLINDER AND DIVERGENCES 



We have considered cylinders of polygonal shapes in this paper; we believe that we have 
covered all cases in which the eigenvalues can be explicitly given. 1 Furthermore, we have 
given numerical results for Dirichlet boundary conditions for a variety of right triangles. It 
would, of course, be of great interest to consider general polygons, in particular regular poly- 
gons, but that is not possible analytically with the exception of the equilateral triangle and 
the square. (For example, only some of the eigenmodes of a hexagon are known.) Instead, 
one may use the method described in Sec. IVIIII to obtain good numerical approximations 
for polygons of arbitrary shape. 



Of course, there are well-known analytic results for circular cylinders [10|, lUl |27|, 40 ]. 
In these cases, however, both interior and exterior contributions to the cylindrical shell are 
included, which is necessary, because otherwise a finite energy cannot be calculated. Put 
another way, the a 2 heat kernel coefficient proportional to the cube of the curvature cancels 



only when interior and exterior modes are included [4l|, |42j. If a 2 7^ 0, the Casimir energy 
cannot be unambiguously computed, because an arbitrary logarithmic scale will add to it. 

Why does this problem not arise here, where it is impossible to compute the exterior 
modes? Let us consider the heat kernel for the square, where the mode sum is especially 
simple. The heat kernel, which lies at the root of the analytic approach followed here, is 
expressed in terms of the operator H = — V 2 

/°° rlh - 1 1 x 

£ J2 e -( fc2+ ^ = ~= J2 ^ lnt - (9-1) 
~°° m,n=l yTTt mn _^ 

The Poisson summation formula (13. 9 p converts the latter sum into a form suitable for small 
t expansion: 

t^'=-W,f^fit^ v "- <<«> 

m=l y y p=l 

The latter sum represents exponentially small corrections as t — > 0+, so by squaring this, 
we get the heat-kernel expansion 

m ~ » - h + i^. < 9 - 3 > 



1 We believe that the analysis in Ref. 29( is in error. 
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where the omitted terms are exponentially small. The three nonzero heat kernel coefficients, 
of course, are proportional to the area, the perimeter, and a constant (this is the generalized 
Weyl theorem): 

ao = d^' ai / 2 = -i£' ai = s^- (9 - 4) 

There is no 02 coefficient. Carrying out the same calculation for the equilateral triangular 
waveguide, we obtain the same result except that a\ = l/(6y/n). These heat kernel co- 



efficients reflect the divergences found in the energy, Eq. (13.141) . as explained in Ref. [42 ]. 
Somewhat mysteriously, the corner divergences, which are physically present, appear in the 
heat-kernel machinery only in the a x coefficient, and not in the a<i coefficient. Because the 
surfaces are flat, there are no curvature divergences as are present in the circular case. In 
general, the corner angles on appear in the coefficient a\ as 

C 



a 1 



48 spii 48 y/n ^ \oti tt 



See Eq. ( 13. 15ft and Refs. [43l-l45| . Thus, the same remarks apply to the other polygonal 



shapes considered here. For a regular iV-polygonal cross section, 



1 JV-1 . , 

ttl = lV^iV^2' (9 ' 6) 

which reduces to the expected result for a circle, 1/(12-*/^), as N — > 00. 

X. PLANE FIGURES 

Perhaps it is even more surprising that finite values are obtained for Casimir energies for 
polygons in a plane, since the corresponding results for a circle are divergent, even when 
interior and exterior contributions are included Q, (47)]. Indeed, Ambj0rn and Wolfram gave 
Casimir energies for rectangles in a plane jjj]. 

The d — version of the Casimir energy ( 12. 2ft is of course 



E 



^J>mn, (10.1) 



which may be immediately evaluated for the four figures considered in this paper via the 
Chowla-Selberg formula, or from the Dirichlet L-series formulae. The results are given in 



Table E (Note there appears to be a transcription error in the formula given in Ref. |28|.) 
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Cross section Dirichlet Neumann A/P 2 



Hemiequilateral Tr. 0.174790/ft -0.238159/fe 0.038675 



Isosceles Right Tr. 0.113080/a -0.202939/a 0.042893 



Equilateral Tr. 



0.0877806//i -0.214519//t 0.048113 



Square 



0.0410406/a -0.220759/a 0.0625 



TABLE I: Casimir energies for plane figures. The second and third columns give the Dirichlet and 
Neumann energies, expressed in terms of the height h for the figures derived from the equilateral 
triangle, and the side a for the figures derived from the square. The fourth column gives the 
dimensionless ratio of the area to the square of the perimeter of the cross section. All results refer 
to interior contributions only. 

These results are plotted versus A/P 2 in Fig. [6j In Fig. [7] we contrast the PFA, given 
here by 



with the Dirichlet Casimir energies found both analytically, and, for general right triangles, 
numerically using the method of Sec. IVIIII 

XI. ANALYSIS AND CONCLUSIONS 

Our results are summarized in Table [TTJ 

We plot the results in Fig. [HJ In both the E and H mode, the results follow a smooth curve 
for all four polygonal cross sections. Even the circle seems to follow the same pattern for the 
E modes, but deviates significantly for H modes, which is hardly surprising since exterior 
physics is included for the circular boundary. We would speculate that if Casimir energies 
for other regular polygonal cross sections could be evaluated numerically, they would lie on 
this universal curve. As we have seen, for the E or Dirichlet modes, the numerical results 
for arbitrary right triangles lie on the same universal curve, which closely approaches that 
expected from the proximity force approximation for small angles. 

Of course, not all cross sections will do so, as illustrated in Fig. [9j This shows the 
trajectory of the Casimir energy for cylinders of rectangular cross section, written as a 
function of the variable A/P 2 . Here, the ratio of the sides of the rectangles can be written 




(10.2) 
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FIG. 6: (Color online) Graph of the interior Casimir energy E for various plane figures, multiplied 
by the square-root of the area A. These are plotted as a function of the dimensionless ratio of the 
area A to the square of the perimeter P, A/P 2 . The points shown are for the 30°-60°-90° triangle, 
the isosceles right triangle, the equilateral triangle, and the square. The upper curve shows the 
result for Dirichlet modes, the lower curve the energy for Neumann modes. 




0.08 p2 



FIG. 7: (Color online) Graph of the interior Dirichlet Casimir energy E for various plane triangles 
and the square, multiplied by the square root of the area A. These are plotted as a function of the 
dimensionless ratio of the area A to the square of the perimeter P, A/P 2 . The points shown are 
for the 30°-60°-90° triangle, the isosceles right triangle, the equilateral triangle, and the square, 
as well as other right triangles evaluated numerically. The curve is the PFA approximation, which 
should become exact when one of the angles of the triangle goes to zero. 
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Cross section Dirichlet Neumann EM A/P 2 

Hemiequilateral Tr. 0.0163745 -0.0204438 -0.00406928 0.038675 
Isosceles Right Tr. 0.0131650 -0.0227063 -0.00954135 0.042893 
Equilateral Tr. 0.0102705 -0.0265477 -0.0162772 0.048113 

Square 0.00483155 -0.0429968 -0.0381653 0.0625 

Circle 0.00193145 -0.0445355 -0.0426041 0.079577 

TABLE II: Casimir energies per unit length for cylinders of various cross sections. The second, 
third, and fourth columns give the E, H, and EM (perfectly conducting) energies/length multiplied 
by the cross sectional area. The fifth column gives the dimensionless ratio of the cross-sectional 
area to the square of the perimeter of the cross section. All results refer to interior contributions 
only, with the exception of the final row, which gives the energies for a shell of circular cross section 
including both interior and exterior modes. 
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FIG. 8: (Color online) Graph of the interior Casimir energy per unit length £ for cylindrical 
waveguides with various cross sections, multiplied by the cross-sectional area A. These are plotted 
as a function of the dimensionless ratio of the area A to the square of the perimeter of the cross 
section P, A/P 2 . The points shown are for the 30°-60°-90° triangle, the isosceles right triangle, 
the equilateral triangle, and the square. The last point is for a circle, including both interior and 
exterior modes. The upper curve shows the energy for E modes (Dirichlet), the lower curve the 
energy for H modes (Neumann), and the intermediate curve the energy for the sum of the two 
modes, that is, for a perfectly conducting cylinder. 
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FIG. 9: (Color online) Graph of the interior Dirichlet Casimir energy per unit length £ for cylinders 
of rectangular cross section, multiplied by the cross-sectional area A. This quantity is plotted as a 
function of the dimensionless ratio of the area A to the square of the perimeter of the cross section 
P, A/P 2 . The upper curve is the interpolation of the triangular and square energies shown in 
Fig. El 



as 



b 
a 



l-8£ 



16 



11.1) 



8£ 2£ 

where £ = A/P 2 , which takes on the maximum value of £ = 1/16 for a square. It will 
be noted that compared to what was found above for the triangles, the behavior of the 
rectangle energies is very different, changing from positive to negative at £ = 0.057902, 
which corresponds to b/a = 1.74437. 
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Appendix A: Chowla-Selberg formula 

Since it is relatively unfamiliar, in this appendix we sketch the derivation of the Chowla- 
Selberg formula from the Abel-Plana formula, 

E /(n) = r dn f(n) + 1/(0) + i [°° d ^'^ , (Al) 
n=o ^ 1 Jo e -1 



which requires that f(z) be analytic in the right half-plane. If f(z) is analytic in the left 
half-plane, we can derive the identical formula for ^^L /(— n). 
Now consider the sum 

oo oo 

s = E E"/( m ^ ( A2 ) 



m=— oo n=— oo 



where the double prime means that the single term m = n = is omitted. We take 

f(m, n) = (am 2 + bmn + cn 2 ) _s , (A3) 
where Re s > 1 so the sum exists. It is easy to show that 

oo oo oo oo 

S = 2 E E n ) + 2 E E -™) - 2a " S C(2s) + 2c- s C(2s). (A4) 

m=l n=0 m=l n=0 

We apply the Abel-Plana formula to the two sums in Eq. (IA4I) . with the result 

S = 2cT s C(2s) + II + III, (A5) 

where 

°° poo 

II = 2E / dn [f(m, n) + f(m, -n)} , (A6a) 

m=l J ° 

and 

/■OO 1 00 

III = 2i / eft — E **) " - lt ) + Z( m ' - ' ( A6b ) 

■'o e - 1 m=1 

where the different terms have branch points in the 2nd, 3rd, 1st, and 4th quadrants, re- 
spectively. 

It is rather easy to evaluate the integral II in terms of the integral 



/■oo 

/ dt(t 2 + 1) 
Jo 



v^r(s - 1/2) 



2r(a) 

The end-point contributions cancel, with the result 



(A7) 



11 = v^^CPs - (A8) 
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where A = 4ac — b 2 . 

The final integral III is evaluated by shifting the contour by an imaginary amount (the 
shifts along the imaginary axis vanish) and encircling the branch point, so that we obtain 

trr\ 1_2s 00 00 poo 

— ] eT s ^ m x ~ 2s J^(-2) / du{e 1Tis -e- ivs ){u 2 -l)- s e-^ up/c 
° / m=l k=l 



fg—imkbTr/c _|_ ^imkbn/c^j 
is/2 1 

r c T(s) 



ai/4-s/2 1 00 

2 s + s/2^ — __ cos(7r^/cy- 1 /V 1 _ 2s (p)/C-i/ 2 ( 7 rpv / A/c), (A9) 

p=l 



which uses the representation of the modified Bessel function, 



Kv (z) = ^ {Z,2) " / dt (t 2 - iy- l ' 2 e- zt . (AlO) 



T{v + l/2) J, 

Using Eqs. (TA~8l) and fTA^j) in Eq. fTAol) . we find the Chowla-Selberg formula (13. 2p with a <-»■ c 
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